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Abstract. Recent evidence has emerged that the Cepheid PL relation in the LMC is non- 
linear in the sense that the existing data are more consistent with two lines of differing slope 
with a break at a period of 10 days. We review the statistical evidence for this, the impli- 
cations for the extra-galactic distance scale and CMB independent estimations of Hubble's 
constant and briefly outline one possible physical mechanism which could cause this non- 
linearity. 
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1. Introduction 

The vast majority of work on the Cepheid 
based extra-galactic distance scale relies on the 
assumption of a linear PL relation in the LMC 
in both y and / bands: Myj - a + blogP. 
Recent evidence has emerged, based on a va- 
riety of Cepheid data, that in the LMC, the 
Cepheid data are more consistent with two 
lines of differing slopes separated at a period 
of 10 days. CuiTent data indicates that this non- 
linearity is present in the B,V,I,J,H bands and 
is marginally linear in the K. Here by PL rela- 
tion, we mean the Cepheid PL relation at mean 
light. 

However, we can also investigate the PL re- 
lation at different phases by phasing all light 
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curves to a common starting epoch. Figure 1 
displays multiphase PC (left panel) and PL 
(right panel) relations in V - / colors and 
V & I bands, respectively, evaluated in this 
way using OGLE data for LMC Cepheids. The 
left panel of figure 1 also presents Galactic 
Cepheid PC relations. In order to calculate 
these, all observations are evaluated relative 
to a common starting epoch. These data have 
been corrected for reddening exactly as sug- 
gested by the OGLE team. The dynamic nature 
of both PC/PL relations are evident as is the 
fact that changes in one are usually reflected in 
the other At a phase of 0.82, there is a dramatic 
break in both V and / band PL relations: per- 
haps the strongest evidence yet of nonlinearity 
in the Cepheid PL relation. Phase 0.0 is at max- 
imum V band light. We see that at this phase 
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the PC relation is very flat. This is similar to 
Galactic PC relations at maximum light. The 
mean light relation is the average of such rela- 
tions over phase and so it may be expected that 
the behavior at phase 0.82 has some effect on 
the mean light PL relation. Our approach thus 
reUes on studying PC/PL relations as a func- 
tion of phase. 

A nonlinearity in the form of a quadratic 
Cepheid PL relation has been suggested be- 
fore (Marconi et al., (2005) and references 
therein). However our approach is different in 
that the multiphase plots suggest two straight 
lines joined at a period of 10 days and the phys- 
ical mechanism we invoke to bring this about 
has not been studied before in the literature. 

2. Statistical tests and the 
nonlinearity of the PL relation 

In order to test for nonlinearity, statistical tests 
are required. Figure 2 displays two artificial 
PL relations, with and without a break at 10 
days. In both cases, the intrinsic scatter around 
the basic regression line is added to mimic 
the real width of the instability strip. By eye 
it is very difiicult to determine nonlinearity, 
yet an F test applied to these data correctly 
determines the intrinsically nonlinear PL re- 
lation. Further it is insuflicient to simply test 
for linearity/nonlinearity by simply comparing 
the short/long period slope and their associated 
standard deviations. Estimations of slope and 
their standard deviations such as "the slope is 
X ± 6x" implies the probabiHty that the slope 
is in the interval {x - 6x,x + 6x) is 1 - or, 
where a is the desired significance level. If 
A is the event that the short period slope is 
wrong and B is the event that the long period 
slope is wrong, P(A) = a and P(B) = a. Then 
in comparing the long and short period slopes 
purely by a comparison of the standard devia- 
tions, the probability of at least one mistake is 
PiA [JB) = PiA) + P{B) - P{A r\B) = 2a-a^. 
If 1 > or > 0, then 2a - aP' > a which is the 
probability the F test makes a mistake. 

Kanbur and Ngeow (2004), Ngeow et al 
(2005), Kanbur and Ngeow (2006), Kanbur et 
al (2007), Koen et al (2007) appUed the F 
test, testimator, Bayesian information methods. 



non-parametric procedures amongst others, to 
the OGLE LMC Cepheid data, MACHO LMC 
Cepheid data, 2MASS IR data cross-correlated 
with MACHO, OGLE LMC data appended 
with data from Sebo et al (2002), Caldwell and 
Laney (1991) and Persson et al (2004) to study 
the linearity/nonUnearity of the LMC Cepheid 
PL relation at mean light at wavelengths rang- 
ing from B to K. All these data sets and statis- 
tical procedures yielded the result that current 
data suggest the LMC Cepheid PL relation is 
nonlinear at a period of 10 days from B to H 
whilst the K band PL relation is marginally lin- 
ear 

One possible objection to this is that there 
are reddening errors particularly for the longer 
period LMC Cepheids ie. the reddening error 
is a function of period. A reddening error with 
period sufficient to make the mean light LMC 
PL relation Unear in an F test will result in 
an extinction error as a function of period: if 
this same extinction error is applied to maxi- 
mum light, it will imply that LMC Cepheids 
get hotter/bluer at maximum light as the pe- 
riod increases. Such behavior is in stark con- 
trast to Galactic Cepheids which display a flat 
PC relation at maximum light. There is credi- 
ble physics explaining the behavior of Cepheid 
PC relations at maximum Ught: the interaction 
of the stellar photosphere and hydrogen ion- 
ization front (HIE, Kanbur Ngeow and Buchler 
2004). 

Another possible objection is that there is a 
lack of long period data. While the number of 
short period Cepheids is indeed much greater 
than the number of long period stars, the F test 
takes into account both the number and nature 
of the data. Thus datasets with fewer numbers 
of long period Cepheids determine the long pe- 
riod slope with a higher error which in turn 
makes it harder for the F test to be significant. 

3. Effect on the distance scale and Ho 

Ngeow and Kanbur (2006c) calibrate the 
Hubble diagram of SNla supernovae and es- 
timate H() using linear and nonlinear LMC 
Cepheid PL relations as calibrators. If the true 
relation is indeed nonlinear and a linear rela- 
tion is used, an error of 1 - 2% in //q results. 
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Fig. 1. Left panel: Multiphase PC relations in V-I colors for extinction corrected OGLE LMC and Galactic 
Cepheid data. Right panel: Multiphase PL relations in V, I for extinction corrected OGLE LMC Cepheid 
data. These figures are adopted from Ngeow and Kanbur {2006a). 



This may seem small but with significant ef- 
fort being made to reduce zero-point uncertain- 
ties (Macri et al., (2006)) and important bene- 
fits accruing from an independent Cepheid dis- 
tance scale accurate to less than 5%, an un- 
derstanding of this possible PL nonlinearity is 
important. Moreover if this nonlinearity is in- 
deed real, then it will be vital to understand 
its physical cause for theories of stellar pulsa- 
tion/evolution. 

4. A Possible Physical Mechanism to 
explain this nonlinearity 

We postulate that the same physical mecha- 
nism responsible for the flat/flatter PC relations 



at maximum light for Galactic and Magellanic 
Cloud Cepheids affects the properties of the 
mean light PC relation. The HIF and stellar 
photosphere move in the mass distribution as 
the star pulsates but they are not comoving. 
The distance between them, in terms of the 
mass distribution, can vary according to phase, 
period and metallicity. If the hydrogen ioniza- 
tion front (HIF) and stellar photosphere are in 
contact at low densities, then the temperature 
of the photosphere is related to the temperature 
at which hydrogen ionizes. However the prop- 
erties of the Saha ionization equation relevant 
to Cepheid envelopes means that the tempera- 
ture at which hydrogen ionizes is relatively in- 
sensitive to temperature at low densities. This 
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Fig. 2. Simulated PL relations, one of which is intrinsically linear and the other one which is intrinsically 
nonlinear. They look similar and can only be distinguished by statistical tests. This figure is adopted from 
Ngeow and Kanbur (2006b). 



can change the PC properties of the Cepheid 
envelope. The amount of this change is related 
to period, phase and metalUcity. Since changes 
in the PC relation affect the PL relation, we 
also postulate that this mechanism can explain 
the nonlinear LMC PL relation: preliminary 
full ampUtude model calculations support this 
idea. 



period for a sample of RRab stars using data 
taken by Cacciari et al (2005). Open/closed cir- 
cles/triangles overtone/fundamental pulsators 
for M3/M15 respectively. For the fundamental 
mode pulsators of M3 we clearly see a flat PC 
relation at minimum light. This flatness occurs 
due to the nature of the HIF-stellar photosphere 
interaction. 



5. Extension to RR Lyraes 

The same physical mechanism can also effect 
the PC relation of RR Lyraes and is in fact 
responsible for the flat PC relation found for 
RRab stars at minimum light in a number of 
colors. The different LjM ratios and hotter ef- 
fective temperatures change the nature of the 
HIF-stellar photosphere somewhat but the es- 
sential idea is the same. 

Figure 3 presents a plot of (B - V) and 
(y - /) color at minimum light against log 



6. Conclusions 

We review recent evidence for a break in the 
LMC Cepheid PC/PL relation. We conclude 
that statistical tests are important and that cur- 
rent data indicate strongly that a nonUnearity is 
present in the LMC Cepheid data at a period of 
10 days. The multiphase relations and consid- 
eration of PC relations at maximum light im- 
ply that it is unlikely that these nonlinearities 
are due to reddening errors. However it will be 
important to further test a more uniform sam- 
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Fig. 3. RR L5Tae colors at minimum light as a function of period for M3 and M15. Overtone/fundamental 
pulsators are open and closed circles/triangles . 



pie of LMC Cepheid data, for example, from 
the LMC Shallow Survey (Gieren 2007). We 
postulate that this break in the Cepheid PC and 
hence the PL relation is caused by the inter- 
action of the HIF and stellar photosphere. The 
variation of this interaction with period, pulsa- 
tion phase and metalUcity may lead to the ob- 
served characteristics of the PL relation in the 
Galaxy and Magellanic Clouds. We also note 
that the same physics can explain the behavior 
or RR Lyrae PC relations at minimum light. 
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